Forming Nonsingular Black Holes from Dust Collapse 
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The dynamics of the gravitational collapse is examined in the realm of string based formalism 
of D-branes that encompass General Relativity as a low energy limit. A complete analytical so- 
lution is given to the spherically symmetric collapse of a pure dust star, including its matching 
with a corrected Schwarzschild exterior spacetime. The collapse forms a black hole (an exterior 
event horizon) enclosing not a singularity but perpetually bouncing matter in the infinite chain of 
spacetime maximal analytical extensions inside the outer event horizon. This chain of analytical 
extensions has a structure analogous to that of the Reissner-Nordstrom solution, except that the 
timelike singularities are avoided by bouncing barriers. The interior trapped bouncing matter has 
the possibility of being expelled by disruptive nonlinear resonance mechanisms. 

PACS numbers: 



Introduction - Black holes are solutions of vacuum Gen- 
eral Relativity equations describing the exterior space- 
time of the final stage of gravitationally bounded sys- 
tems whose masses exceeded the limits for a finite equi- 
librium configuration Geometrically a black hole may 
be described as a region of asymptotically flat space- 
times bounded by an event horizon hiding a singularity 
formed in the collapse. Fundamental theorems by Is- 
rael and Carter state that the final stage of a general 
collapse of uncharged matter is typically a Kerr black 
hole, which has an involved singularity structure. Nev- 
ertheless, for a realistic gravitational collapse we have 
no evidence that the Kerr solution describes accurately 
the interior geometry of the black hole. On the contrary, 
the best theoretical evidence presently available indicates 
that the interior of the black hole thus formed is analo- 
gous to the interior of a Schwarzschild black hole with a 
global spacelike singularity @. The simplest way of form- 
ing such structure is by the spherical collapse of dust, as 
originally shown in the classical paper of Oppenheimer 
and Snyder However, as singularities cannot be em- 
pirically conceived, this turns out to be a huge pathology 
of the theory. 

Albeit the cosmic censorship conjecture there is no 
doubt that the General Theory of Relativity must be 
properly corrected or even replaced by a completely new 
theory, let us say a quantum theory of gravity. This de- 
mand is in order to solve the issue of the presence of sin- 
gularities predicted by classical General Relativity, either 
in the formation of a black hole or in the beginning of the 
universe. While a full quantum gravity theory remains 
presently an elusive theoretical problem, quantum grav- 
ity corrections near singularities have been the object of 
much recent research, from loop quantum cosmology^ 
to string-based formalism of D-branes theory Q- In the 
latter scenario extra dimensions are introduced consti- 
tuting the bulk space. All matter would be trapped on a 
4-dim world-brane spacetime embedded in the bulk and 



only gravitons would be allowed to move in the full bulk. 
At low energies General Relativity is recovered || but at 
high energy scales significant changes are introduced into 
the gravitational dynamics and the singularity can be re- 
moved. 

In this context, we examine the gravitational collapse 
of a pure dust star in the realm of braneworld theory 
in a 5-dim bulk, a procedure analogous to that of Op- 
penheimer and Snyder Q. We reproduce the interior 
Friedmann- Robertson- Walker solution and obtain its an- 
alytical matching with a corrected Schwarzschild exterior 
geometry. The maximal analytical extension of such ge- 
ometry is analogous to that of the Reissner Nordstrom 
solution but now with no charge and no singularity at 
all. Our configuration exhibits at most two event hori- 
zons enclosing perpetually bouncing matter with no sin- 
gularity formation. We also examine the implications 
of such corrections in the context of the General Rel- 
ativity tests, and in the Hawking evaporation of black 
holes. We rely basically on Refs. @ for the modified 
Einstein's equations in a 4-dim braneworld embedded in 
a 5-dim bulk. Related literature on the problem of non- 
singular black holes is properly addressed at the end of 
the paper. Throughout the paper we use units such that 
Gjv = c = 1; however for clarity we maintain Gm in all 
expressions. 

The Model - We assume a spacetime braneworld model 
embedded in a 5-dim bulk with a timelike extra dimen- 
sion, whose matter content is a spherically symmetric 
collapsing dust with density p. In a coordinate system 
comoving with dust the interior geometry is still shown 
to be a Friedmann-Robertson- Walker metric fioj 
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with its dynamics given by the first order modified Fried- 
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maim equation 
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where Eo is a constant of motion associated with the 
dust density, Eq = pa 3 , and a is the negative brane 
tension (TTj. Gat is Newton's constant on the brane. As- 
suming initial conditions for the collapse d(0) = and 
a(0) = 1, we get 
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The last term in the RHS of Equation is a correction 
due to the bulk-brane interaction that eventually avoids 
the singularity, acting as a bouncing barrier. Its negative 
sign is consequence of the timelike character of the extra 
dimension. The scale of the bounce is given by p ~ 2\a\ 
that corresponds to a minimum of the scale factor a m ~ 

/ \l/3 

(Eo/2\a\) . 

We determine the spherically symmetric metric outside 
the collapsing star from its matching to the metric |T]), 
at the surface defined by r = 7 in comoving coordinates. 
To this end let us transform the comoving coordinates of 
([T]) to Schwarzschild coordinates (T, R, 9, 4>) through the 
equations 



T = F[S(r, £)], R=a(t)r, 
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with F[S(r,t)] an arbitrary function of S, where S(r,t) 
satisfies 
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The line element (JTJ is diagonal in this new coordinate 
system and we automatically guarantee that 
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where the second equality follows from Equation (J2J). It 
is straightforward to verify that the solutions for are 



S(r,t) =S + a(l 
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and 



S(r,t) = S + C 
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where S, a and C are arbitrary constants. With the 
use of ([2]) the above integrals arc analytically soluble for 
k > 0. For the physical domain of parameters to be 
considered here, S turns out a monotonous function of a 
which can be properly inverted to express a in terms of 
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FIG. 1: Plot of the polynomial V{R) = R grr for dust masses 
M < M» (no black hole), M = M» (extremal black hole) 
and M > M, (black hole with outer horizon R+ and inner 
horizon _R_). The figure corresponds to \a\ — 0.05, in units 
G N = c = 1. 



5, for an explicit range of a. The remaining task is to 
choose F in terms of iS. A choice can be suitably made 
for the case k > [13] so that at the surface of matching 
r = 7 we obtain 
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is the total mass of the collapsing dust. The equation of 
the surface of the star is given by R = a(t) 7. In the 
remaining of the paper we are restricted to k > 0. In 
this instance the motion of dust is oscillatory, bouncing 



between a = 1 and 



the two real roots of a = 



Eq. ©. The assumed initial conditions for the collapse, 
a(0) = 1 and d(0) = 0, require then that a m < 1, leading 
to a restriction on a, namely \a\ > 2Eq- 

Maximal Extension - In order to examine the maximal 
extension of the exterior geometry @, we need to know 
whether, and under what circumstances, the configura- 
tion forms event horizons. A necessary condition is that 
the mass M of the collapsing star equals or exceeds a 
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(\a\ fixed), namely, that the polynomial V{R) = R A grT 
has one i?*, or two < R + ) roots respectively (cf. Fig. 
1). Otherwise we cannot have formation of event hori- 
zons. The maximal analytical extension for the space- 
time (9) is considered for condition M > M* (cf. (11)) 
and 7a m < R- . It is analogous to that of a Reissner- 
Nordstrom black hole[lH for Q 2 < M 2 but without 
timelike singularities (cf. Fig. [2] and its caption). It 
is worth remarking that the surface of the collapsing 
dust - once crossing R + - must necessarily cross i?_ , 
namely, ya rn < R_ , so that a stable black hole forms 
with trapped perpetually bouncing matter. We should 
mention that extensive computer tests of the matching 
(|9j) with parameters taken in a physically satisfactory 
domain (basically in accordance with (|11|) . and satisfying 
the lower bound condition for \a\ in order to have oscilla- 
tory motion) necessarily lead to 7a m < i?_ , assuring that 
Eq. (jlip is in fact a consistency condition for the forma- 
tion of a nonsingular black hole with two horizons. For 
illustration let us consider S = 1 and a — 0.01 in ([7]). Wc 
also take E = 1.04 x 10~ 21 mT 2 (corresponding to one 
thousand times the solar mass density), 7 = 6.96 x 10 8 m 
(corresponding to one solar radius) and a = 10 s (a large 
value to tentatively approach General Relativity). We 
obtain for S(a) a strictly monotonous function in the 
non-negative range of a, and i?_ > ja m - specifically 
ja m = 0.12058048347 and i?_ = 0.12058048511. As \a\ 
decreases the difference (R-— r ya m ) increases by the same 
order of magnitude. 

Summarizing, the gravitational collapse of a dust star 
with mass greater than M* produces a black hole - de- 
fined as the region of spacetime enclosed by the outer 
event horizon R + — enclosing not a singularity but per- 
petually bouncing matter in the infinite chain of space- 
time maximal extensions inside the outer event horizon. 

We note that the critical mass depends solely on the 
parameter a. To have an idea of the order of mag- 
nitude of a for event horizon formation, let us take 
M* ~ IAMq, the Chandrasekhar limit. This yields 

cr*| ~ 10 -2 kmT 2 . A star with the Chandrasekhar mass 
would not form an event horizon if the brane tension is 
smaller than 10 -2 km~ 2 , so that this value establishes 
a lower bound for \a\. It is also worth noticing that the 
geometry ^ is a solution of brane corrected Einstein's 
equations R a b + E a (, = 0, where E a b is the traceless pro- 
jection of the Weyl curvature on the brane[9j; we can 
easily extend Birkhoff theorem for these equations. 

Observational Consequences and Final Comments - 
For the black hole geometry ((9]) we have derived the two 
classical tests of General Relativity, the planetary perihe- 
lia precession and the bending of light. In a perturbative 
procedure [l4T| we obtain the perihelion advance per rev- 
olution 

6irG 2 N M 2 / _ 3G N M 2 _ 15G^M 4 \ 
h 2 V 4tt|ct|/i 4 ~ 8tt|ct|/i 6 /' 

with h 2 — GffAla(l — e 2 ), where e is the eccentricity of 
the orbit and a its semimajor axis. The first term on the 




FIG. 2: Penrose diagram for the spacetime Q assuming M > 
M* (cf. |H])). The infinite chain of asymptotically flat regions 
I (00 > R > R+) are connected to regions III (R- > R > 
7a m ) by regions II (7?+ > R > i?_). The dashed portion 
of regions III bounded by R = *ya m and R = is not part 
of the spacetime. Once crossing R- (the Cauchy horizon) 
the matter remains perpetually bouncing between R = 7 and 
R ~ ya m . Dashed arrowed lines illustrate timelike geodesies 
of dust particles and observers bouncing at R = 7<i m . 



RHS corresponds to the General Relativity correction. 
Analogously for the bending of light rays passing in the 
neighborhood of a spherical massive body, the deflection 
angle of the asymptotes is given by 

_ AG N M / _^5_ M \ 
R V 256 |a|i? 3 / 

where R is the radius of the body. Again the first term 
on the RHS corresponds to the General Relativity cor- 
rection. For M — Mq the corrections to the effects 

predicted by General Relativity are neglegible for \a\ 
above the lower limits previously evaluated. However 
these corrections might be important for masses and/or 
scales much larger than that of the solar system. 

In the case of Hawking radiation [lEj], the calculations 
give qualitatively analogous results, with a Planckian 
thermal spectrum of created particles, but with a cor- 
rected Hawking temperature depending on M and |er|, 

(R+ - R-) 

T » ^ ^rTIT- < 12) 
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We note that in the extremal case we have R + = i?_ = 
3GnM/2 implying that Th — * continuously as — » 
J2_. The constants A and B have the expression A = 
(R+ + R-) - 2G N M and B = 3G N M 2 /{Air\a\R + R-). 
The observation of Hawking radiation could, in principle, 
allows to test our results for finite \cr\. 

Black Hole Thermodynamics - In this section, some 
thermodynamical results are derived for the case of quasi- 
extremal black holes. In this case the roots of our poly- 
nomial P(R) are given by R± ~ 3GnM/2 ± e, where e is 
a small deviation from the extremal case and given by 
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By defining the outer horizon area as A onteT 
obtain that 
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calculated about the extremal configuration(M», |cr* |). In 
this approximation the Hawking temperature (|12p as- 
sumes the form 



T, 
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the substitution of which in (fl"3l) results 
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We can therefore associate the horizon area of the quasi- 
extremal black hole with the entropy S = Amter/4Gjv, in 
accordance to Bekenstein's definition fl7|, turning Equa- 
tion Q14[) into an extended Second Law of Thermodynam- 
ics, with an extra work term connected to the variation 
of the brane tension. 

Final Comments - Finally, according to previous work 
on perturbations in bouncing matter [lfjj, nonlinear reso- 
nance mechanisms typically may occur in the collapsed 



bouncing matter due to its parametric fluctuations as 
produced for instance by a perturbative scalar field. Such 
a mechanism can turn the black hole configuration un- 
stable leading to a disruptive ejection of mass. In fact, 
the dynamics of oscillatory perturbations of the interior 
spacetime geometry and matter content present small do- 
mains in the parameter space where the dynamics is res- 
onant. The remaining parametric space corresponds to 
stable configurations. The parameters \a\ and the am- 
plitude of the perturbations regulate the driving of the 
system from one region to the other. We will approach 
this issue in a future publication. 



Also the presence of Cauchy horizons R- in the maxi- 
mal analytical extension of the geometry §§§ poses a ques- 
tion on the stability of the spacetime. Due to its simi- 
larity with the Reissner-Nordstrom spacetime we should 
expect that the flux of energy of test fields may diverge 
on crossing i?_ (cf. fl3j and references therein). How- 
ever a complete treatment of the problem should include 
higher-order nonlinear terms to provide sufficient condi- 
tions for instability. This work is in progress and will 
also be the object of a future publication. 



An analogous problem - the end state of a collaps- 
ing homogeneous scalar field - was attacked by Bojowald 
et al. [18( in the realm of loop quantum cosmology. 
In their paper the singularity avoidance basically arises 
from semiclassical corrections in the geometric density 
I /a 3 , which does not diverge as a — > 0. However the 
interior spacetime containing the scalar field distribution 
cannot be directly matched to a static Schwarzschild ex- 
terior, leaving therefore the dynamics of horizon forma- 
tion and its structure undetermined. A critical threshold 
scale for horizon formation is obtained, indicating the 
possibility of formation of dynamical horizons, but the 
authors were not able to compute the associated criti- 
cal mass because the exterior dynamics remains undeter- 
mined. Their claim that both horizons will become time- 
like and evaporate is not sustained by the calculations in 
our model. Also several proposals for the construction of 
nonsingular black holes in General Relativity appeared 
in the literature. However none of them are solutions 
of Einstein's equations neither are they associated with 
known physical sources (cf. [l9j and references therein). 
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